Synteza termojądrowa

Współczesny model gwiazdy domyka obliczenie tempa produkcji energii
w reakcjach syntezy termojądrowej i powiązanej z nimi produkcji neutrin.

co do zasady wzór  dobrze wyjaśnia źródło energii

cztery atomy wodoru przekształcają się w atom helu

masa atomu helu/cząstki  jest mniejsza niż masa 4 atomów wodoru/protonów

różnica masy  przekształcana jest na fotony 
i neutrina elektronowe

neutrina z prędkością światła uciekają od razu, dlatego odejmuje się je od tempa produkcji energii (dla Słońca jest to tylko 2%, ale dla presupernowej praktycznie 100%)

p=IsotopeData["Properties"]

u=QuantityMagnitude[Quantity[1,"AtomicMassUnit"]*Quantity[1,"SpeedOfLight"]^2,"MeV"]

nuc=IsotopeData[];

num=Dimensions[nuc][[1]]

IsotopeData["He4","AtomicNumber"]

IsotopeData["He4","Stable"]

{IsotopeData[nuc,"AtomicNumber"],IsotopeData[nuc,"BindingEnergy"]}//Transpose; (* 1000x faster than Table below *)

(*Table[{IsotopeData[nuc[[k]],"AtomicNumber"],IsotopeData[nuc[[k]],"BindingEnergy"]},{k,1,num,1}];*)

d={IsotopeData[nuc,"AtomicNumber"],IsotopeData[nuc,"NeutronNumber"],-IsotopeData[nuc,"BindingEnergy"],IsotopeData[nuc,"Stable"]}//Transpose;

d[[-1]]

stable=Select[d,#[[-1]]==True&][[All,{1,2,3}]];

Length[stable]

stable[[1;;12]]

stableNeg=MapAt[Minus,stable,{All,3}];

BE1=ListPointPlot3D[stableNeg,Filling->Bottom,PlotRange->{{0,All},{0,All},Automatic},AxesLabel->{"Z","N","E [MeV]"},
Ticks->{Automatic,Automatic,{7.5,8,8.5,8.8}},BoxStyle->Thick]

BE2=ListPointPlot3D[stableNeg,Filling->Bottom,PlotRange->{{0,All},{0,All},{7.5,8.8}},AxesLabel->{"Z","N","E [MeV]"},
Ticks->{Automatic,Automatic,{7.5,8,8.5,8.8}},BoxStyle->Thick]

BE3=ListPointPlot3D[stable,Filling->Top,PlotRange->{{0,All},{0,All},{-8.8,-7.5}},AxesLabel->{"Z","N","E [MeV]"},
Ticks->{Automatic,Automatic,{7.5,8,8.5,8.8}},BoxStyle->Thick]